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Abstract 

In the framework of a newly developed factorization formalism which is based on NRQCD, explicit 
cancellations are shown for the infrared divergences that appeared in the previously calculated hadronic 
annihilation decay rates of P-wave and D-wave heavy quarkonia. We extend them to a more general case 
that to leading order in v 2 and next-to- leading order in a s , the infrared divergences in the annihilation 
amplitudes of color-singlet QQ( 2S+1 Lj) pair can be removed by including the contributions of color- 
octet operators QQ( 2S+1 (L - l) r ), QQ( 2S+1 (L - 3).///), • ■ ■ in NRQCD. We also give the decay widths 
of Dj — > LH at leading order in a s . 



Since the discovery of charmonium and bottonium, heavy quarkonium physics has drawn much 
attention from theorists. Particularly, studying the decay rates of heavy quarkonium states into light hadrons 
may provide very interesting tests of perturbative quantum chromodynamics (QCD). In the early "color 
singlet model" analysis, quarkonium is considered as a bound state of color singlet QQ pair which is in a 
fixed angular-momentum state. And in the formula of the decay rates, all long-distance nonperturbative part 
was assumed to be factored into the nonrelativistic wavefunction of color singlet QQ or its derivative at the 
origin, and the perturbative part is related to the annihilation of color-singlet QQ which can be calculated 
using perturbative QCD. In the nonrelativistic limit, this early factorization formalism was supported by 
explicit calculations for S-wave decay at next-to- leading order in a s Q . But in the cases of P-wave @ and 
D-wave Q quarkonium decays, the infrared divergences occur, which are associated with the logarithms of 
binding energy in the perturbative calculations of color-singlet QQ annihilation amplitudes. These indicate 
that the decay rates are sensitive to nonperturbative effects beyond those related to the wavefunction of 
color-singlet QQ pair or its derivative at the origin, and not all nonperturbative effects can be factored into 
the color-singlet component of quarkonium. The early factorization formalism is incomplete and fails to give 
correct results in some cases. 

Recently, Bodwin, Braaten and Lepage have developed a rigorous factorization formalism [|J which 
is based on an effective field theory, nonrelativistic QCD (NRQCD). This factorization formalism gives a 
systematic analysis for the decay and production of heavy quarkonium to any given order in v 2 and a s . In 
this factorization formalism, the decay rates can be written as a sum of a set of long-distance nonperturbative 
matrix elements, each of which multiplied by a short-distance perturbative coefficient, which can be calculated 
in perturbative QCD. The nonperturbative factors in this formalism are defined rigorously in NRQCD. This 
factorization formula applies to the decay of S-wave, P-wave equally as well as to the decay of any high 
orbital angular-momentum states. The essential point of this formalism is that quarkonium is treated as a 
state consisting of QQ, QQg and other high Fock components rather than only consisting of color-singlet 
QQ component. In some cases the color-octet component even give dominant contribution to the decay. 

In NRQCD the annihilation rate of a quarkonium state H to light hadrons {LH) can be written as 



where the sum is over all possible local 4-fermion operators O n that annihilate and create a QQ pair, and 
d n is the scaling dimension of O n . 

Since all nonperturbative effects are factored into the matrix elements, the coefficients must be 
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infrared finite and can be calculated in perturbative QCD. The coefficients associated with S-wave decay 
to next- to-leading order both in a s and v 2 have been given in Ref. |4|, [|. There the calculation of QQ 
annihilation rates in full perturbative QCD does not involve infrared divergences, and in deriving the color- 
singlet coefficients one does not need to consider the contribution from color-octet operators in NRQCD, 
which is suppressed by w 4 as compared with that from dominant color-singlet operator in the S-wave decay. 
But for the P-wave and D-wave decays, the color-octet component may contribute at the same order in 
v 2 as the color-singlet component, and we would get infrared divergence coefficients if only color-singlet's 
contribution is considered in NRQCD. ft is very interesting to know how the cancellation of these divergences 
occurs by introducing the color octet operators in NRQCD and get infrared finite coefficients. In this paper 
we will show the explicit cancellation of infrared divergences by analysing the new factorization formulism 
to P-wave and D-wave decay processes. The decay rates of D-wave triplet states will also be given to leading 
order in a s . 

It is known that the coefficients Imf n (a s ) in (|l|) can be determinated by matching the imaginary part 
of QQ forward scattering amplitude in the full theory with that in NRQCD, where QQ is in the corresponding 
color and angular momentum state. In the full theory we recalculate the scattering amplitudes of color-singlet 
P-wave and D-wave QQ in dimensional regularization scheme and get 
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Here for convenience we only give their infrared divergence part and C( 2S+1 Lj) represent finite terms. We 
find when making a substitution -j^- — > In^, the above results will coincide with that derived in Q and ||. 

We use dimensional regularization to control infrared divergences and take quarks to be on-shell in order to 
keep gauge invariance. Because conventional NRQCD is treated under the on-shell condition, we must use 
the same regularization scheme in the full theory and NRQCD so that we can get a compatible result. 
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Fig.l Feynman diagram contributing to ImA4( 2S+1 Lj)NRQCD through the color-singlet operator 
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Fig. 2 Feynman diagrams contributing to ImM( 2S+1 Lj)nrqcd through the color-octet operator 
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In the effective NRQCD theory, we can recalculate these annihilation amplitudes by considering the 
contributions of 4-fermion operators. The tree level contributions come from the diagram in Fig.l, where QQ 
annihilate by color-singlet operator. The essential point is that some color-octet operators also contribute 
to the annihilation amplitudes through one-loop diagrams in Fig. 2, which represent the transition between 
color-singlet component and color-octet component by emitting soft gluons. These color-octet operators are 
e> 8 ( 3 <Si), O&CSo), O s { 3 P ), and O s { 3 P 2 ). By an explicit calculation we get 
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where both infrared (IR) and ultraviolet (UV) divergences appear. The UV divergences can be removed 
after the renormalization of color octet operators. IR divergences are due to the nonperturbative effect for 
the transition of QQ by emitting soft gluons. We will sec in the following that these IR diverdences just 
equal to those that appear in the annihilation amplitude calculated in the full theory. At leading order 
the coefficients of color-octet operators can be obtained by simply substituting the color factors for the 
corresponding color-singlet ones in 0] . They read 

(/mA(»S,))„ = (12) 

(/m/sCS.)). = " {N {-* W - , d3) 
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Substituting the leading order coefficients (|12-15) into (|7|-|l 1|) and comparing them with the results in the 
full theory we find that the infrared divergences are cancelled and finite coefficients for color-singlet 

operators of P-wave and D-wave decays are obtained. It is interesting to note that IR divergences in the 
annihilation amplitudes of color-singlet P-wave QQ can be factored into the matrix elements of color-octet 
S-wave operators, and in the D-wave cases IR divergences can be factored into the matrix elements of color- 
octet P-wave operators, and at order of a% the color-octet S-wave operators do not contribute to the IR 
divergence parts of annihilation of D-wave QQ pair. 
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Fig. 3 Examples of real-gluon emission in quarkonium decay at order o^. Here I represents the 

4-momentum of a soft gluon. 

We can extend this conclusion to a more general case that at leading nontrivial order in v and at 
order af, the IR divergences appearing in the color-singlet QQ( 2 +1 Lj) annihilation amplitude can just be 
cancelled by the contribution of color-octet operators 0( 2S+1 (L — 1),/'), 0( 2S+1 (L — 3)j»), • • •. 

Here we will give a detailed analysis for the cancellation mechanism. We know that at the 
divergences come from only two sources, QQ — > ggg and QQ — > qiqig. We give a specific example of decays 
into 3 gluons (Fig. 3) and the discussion applies also to decays into a light quark-antiquark pair and a gluon. 
In the center of momentum frame, we take the quarks to be on-shell with momentum pi and P2 for Q and 
Q respectively. By energy conservation there must be at least two gluons with momenta of order m, so the 
diagram in Fig. 3(c) contains no infrared divergence. For the diagrams in Fig. 3(a) and (b), the emission 
vertex for the soft gluon with momentum I and one adjacent heavy-quark propagator and one Dirac spinor 
can be approximated as follows: 
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At leading nontrivial order in v 2 we must retain L factors of p in the annihilation amplitudes of the QQ pair 
whose orbital angular momentum is L. The factor of p can come from one of two sources: the purely short- 
distance (infrared-safe) part of the diagram, or the potentially infrared-divergent part ( |l6| ) and (0), which 
consists of the soft gluon and the heavy-quark propagator to which it attaches. If even factors of p come 
from the infrared-divergent part of the diagram i.e., from ( |l6| ) and (0), it is obvious that the divergences 
cancel between the diagrams of Fig. 3(a) and Fig. 3(b). But if odd factors of p from the divergence part, the 
divergences come from the two diagrams add together, rather than cancellation each other. 

In the decay rate, we must integrate the infrared emission factors from the square of the sum of the 
amplitudes over the phase space of gluon. We find results 
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for the case of including one or three factors of p in it respectively. Apart from the infrared emission factors, 
other propagator and vertex factors in Fig. 3(a) and Fig. 3(b) contribute to the short-distance part of the 
decay rate, which is equivalent to the annihilation of QQ through color-octet operators Os( 2S+1 (L — 1),/') 
and 8 ( 2S+1 {L - 3)j») in NRQCD. 

In NRQCD the IR divergence parts of the correspondingly decay rates come from one-loop diagrams 
in Fig. 2 which include color-octet operators 0$( 2S+1 (L — l)jv), 0$( 2S+1 (L — 3)j»), • • •. We can extract the 
IR divergence term through the integration 
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It is obvious that they are the same as PQCD results (|l8|) and Therefore at order of a 3 the source of 

divergence is the transition of color-singlet QQ( 2S+1 Lj) to color-octet QQ( 2S+l (L-l)j>), Os{ 2S+1 (-L— 3)j</), 
• ■ • by emitting one, three • • • soft gluon, and these divergences can be factored into the matrix elements of 
corresponding color-octet operators. 

In order to get a finite result at order a 3 for the decay rates of high angular- momentum quarkonium 
states, we must consider the contribution both from color-singlet and color-octet components. The complete 
formula for decay rate of P-wave states has been given in j(| . In the following we will give the decay widths 
of D-wave triplet only to leading order in a s , and they can be written as 
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where we have neglected the contributions of the matrix elements of color-octet S-wave operators in NRQCD. 
We don't expect the S-wave component to play a major role in the decay od D-wave quarkonium in leading 
order in a s according to the reason in [Q . Here we have used the heavy-quark spin-symmetry, which implies 
that, to leading order in v 2 , the matrix element of 8 ( 3 Pj') in the D-wave triplet can be expressed as one 
nonperturbative parameter 

8 < 1 — |Q 8 ( 3 P 2 )|1— > 

H P - 4 > 

m 4 

which can be determined from lattice calculation or extracted by comparing with experimental data. Here 

we also have used the relation that rate of the matrix element — 1 8 a j is 1 : 9 : 36 for J' = 2 and 

J = 1,2,3; and 1 : 20 for J = 1 and J' = 2, 0. This relation can be derived by considering the rate of El 
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transition from color-singlet 3 Dj to color-octet 3 Pj>, which is proporational to (2 J' + 1) < ^ j j, 



where 
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represents the corresponding 6j coefficient. 



It is interesting to indicate that the ratio for the decay rates of D-wave triplet is R — 76/9 : 1 : 4, 
which is consistent with previous calculations || . The previously calculated R come from the divergence parts 
of decay rate of color-singlet D-wave component, which are now factored into the color-octet P-wave matrix 
elements, and to leading order in a s , we only consider the contributions of color-octet P-wave components. 
Therefore we derive the same R as the previous one. This also happens in the P-wave case Q. 

In the above discussions we have only studied the quarkonium decays at leading nontrivial order 
in v 2 and to next-to- leading order in a s . In principle, it can be applied to any given order in v 2 and a s , 
because it is based on a solid theoretical formalism. But in high order cases the problem is very complicated 
and more operators have to be considered. 
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